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s^ ■ Abstract 

OO 

Q"^ ' This paper studies a class of nonlinear Dirac equations with cubic terms 

Tlj- i in R^^^, which include the equations for the massive Thirring model and 

o 

en 

the massive Gross-Neveu model. Under the assumptions that the initial 
data has small charge, the global existence of the solution in H^ are proved. 
C^ . The proof is given by introducing some Bony functional to get the uniform 

estimates on the nonlinear terms and the uniform bounds on the local 
smooth solution, which enable us to extend the local solution globally in 
time. Then L^-stability estimates for these solutions are also established 
by a Lyapunov functional and the global existence of weak solution in L^ 
is obtained. 



1 Introduction 

We consider the nonlinear Dirac equations 

i{ut + Ux) = -mv + Ni (-u, v) , 

(1.1) 

K'^t -Vx) = -mu + N2{u,v), 
with initial data 

{u,v)\t=o = {uo{x),vo{x)) (1.2) 

where (t, x) G R^, {u, v) G C^, m > 0. The nonlinear terms satisfy the following: 
(Al) A'^i and N2 have following form: 

Ni{u,v) = {aiu + a2U){as\v\ + a^v +a5V), 

N2{u, v) = iPiv + f32v){^3\u\^ + PiU^ + /35m'), 
where a^ and /3fc (A; = 1, ■ ■ • , 5) are constants. 

(A2) For any (-u,!)), 

'Si{mNi{u, v) + ivN2{u, v)) = 0. 

Here and in sequel 3f?Z denotes the real part of Z G C. 

Many physical models like massive Thirrring model and massive Gross-Neveu 
model verify (Al) and (A2). Indeed, for massive Thirrring model ([23]), we have 

(A'"!, A^2) = V(:u,:i;)Q;|Mp|vp = (a-ulf p, af |Mp), a E R^, 



therefore 

U{mNi) = ^{ivN2) = 0. 

For the massive Gross-Neveu model ([17]), we have 

(A^i, A''2) = V(u;v)Oi{uv + uv) = {2av{uv + uv), 2au{uv + uv)), a G R^, 



then 



3fi(mA^i + ivN2) = ^{2ia{uv + uvf} = 0. 



The nonhnear Dirac equation is important in quantum mechanics and general 
relativity ([17] and [23]). There are a number of works on the local and global 
well-posedness of Cauchy problem in different Sobolev spaces, see for examples 
[1], [2], [4], [6], [7], [8], [10], [11], [12], [13], [14], [18], [19], [20], [21], [24] and 
the references therein. The survey of the well-posedness in the nonlinear Dirac 
equation in one dimension is given in [20]. In this paper we consider a class of 
cubic nonlinear Dirac equations in one dimension, which include the equation 
for Massive Gross-Neveu model as an example. To our knowledge, the global 
existence of solution in H^ or in L^ for the Dirac equation of Massive Gross- 
Neveu model is still open [20]. And it was pointed out by Pelinovsky [20] that 
the apriori estimates of the L^— norm (||w(t)||^p -|- ll^^(i^)llLp)"'^''"'' ^^r the equations 
like Gross-Neveu model include nonlinear terms, which may lead to the blow-up 
of the L°° and H^ norms. In addition, the energy of Dirac equation is non- 
positivity . We will overcome the difficulties caused by the nonlinear terms in the 



apriori estimates for p = 2 and the non-positivity of the energy by introducing 
a Bony type functional Q{t) = J J^^y\u(t,x)\'^\v(t,y)\^dxdy, which is similar to 
the Glimm's interaction potential for conservation laws ([16]) and is also similar 
to the Bony functional for Boltzmann equations ([3] and [15]). Then, following 
Bony's approach [3], for any local i7^— solution {u,v) with small initial charge 
/-^(I'^oP + \vo\^)dx, we can get the control on /q /^ \u{x,T)\'^\v{x,T)\'^dxdT by 
using the conservation of the charge J^ooi\'^i^'^)\'^ + \v{t,x)\'^)dx; see Lemma 
2.2. Here the conservation of charge comes from the assumption (A2) on the 
nonlinear terms. These bounds enable us to apply the characteristic method to 
the equations (2.2) to get the uniform L°° bounds on \u\'^ and \v\'^ for extending the 
local //^—solution globally in time; see Lemmas 2.3 and 2.4 for the L°° bounds. 
Then, to study the stability of the if ^—solutions, we have (3.4) and introduce 
Lyapunov functional Li(t) + KQi{t) to derive the L^ estimates on the difference 
between these ii^— solutions. Finally global weak solution in L^ is obtained. We 
remark that Glimm interaction potential was first used by Glimm [16] to establish 
the global existence of the BV solution for the system of conservation laws, then 
was used to study the general conservation laws, and that Bony type functional 
was used to study Boltzmann equations and wave maps in i?^+^, see for instance 
[16], [5], [9], [22], [3], [15]and [25] and the references therein. 
The main results is stated as follows. 

Theorem 1.1 Suppose that {uq,Vo) G H^{R^) and that the norm H^olliacjiji-) + 
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ll^o|li2(j:ji) < Co for some small constant Cq > 0. Then (1.1-1.2) admits a unique 
global solution {u, v) G C([0, cx)); H\R^)) n C\[0, oo); L\R^)). 

Theorem 1.2 Suppose that {uok,vok) E H^{R^) with ||Mofe||i2(Ri) + ||t'ofc||i2(^i) < 
C'o (1 < ^ < oo) such that \\uok — Mooo||l2(_ri) + \\vok — Vooo\\L^{m) — ;■ as k ^ oo 
for some (woooi^ooo) ^ L'^{R^). Let {uk,Vk) be the classical solution of (1.1) 
(1 < k < oo) taking {uok,vok) as its initial data. Then there exists a weak 
solution {uoo,Vryo) G C( [0, +oo), L^(i?^)) of (1.1) such that 

\\Uk — Uoo\\c{[0,T],L'2{R^)) + \\Vk — ^oo||c([0,T],L2(R2)) — )■ 

as k -^ OO for any T > 0. 

2 Global classical solution 

Since C^{R^) is dense in H^{R^), we consider in this section the problem for 
{uq,vq) G C^{R^). (1.1) is a semihnear hyperbolic system. Local existence of 
the solutions in H^ for the semilinear strictly hyperbolic system like (1.1) can 
be proved by the standard methods using the Duhamel's principle and the fixed 
point argument, see for instance [2], [10] and [22]. Now we suppose that for 
{uo,Vo) G C^, the initial value problem (1.1-1.2) has a smooth solution {u,v) G 
C{[0,ny,H\R^)) n C\[(},T,y,L^{R^)) for some T, > 0. Then the support of 
{u{x,t),v{x,t)) is compact for each t G [0,T*). To extend the solution across the 



(2.2) 



time t = T^, we only need to show that 

sup {\Ht)\\H^ + \\v{t)\\H^) < C'(r,)exp(C"(T,)) (2.1) 

0<t<T, 

where C'(A) and C"(A) are polynomials of A independent of T^,. To this end, we 
multiply the first equation of (1.1) by u and the second equation by v, then 

{\u\'^)t + {\u\'^)x = mi{uv - uv) + 2^{iJhu), 
{\vf)t - (kHx = -mi{uv - uv) + 2'R{iN'2v). 
By Assumption (A2), we have 

{\u\^ + \vW + {\u\''-\v\^)^ = Q, 

which leads to the following, 

Lemma 2.1 For any t G [0,T*) there holds 

{\u{t,x)\^ + \v{t,x)\^)dx= / {\uq{x)\'^ + \vQ{x)\^)dx. 



Moreover, by Assumption (Al), we can find a constant c > such that for 
any (-u, v) there hold the following, 

|iVin| + |iV2t;| < c|n|>|l (2.3) 

Then, (2.2) implies that 

{\u\'^)t + {\u\'^)x<rQ{t,x), 

(2.4) 

{\v?)t-{\v?)x<r,{t,x). 
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Here 

ro(t,x) = m(|wp + \v\'^){t,x) + 2c\u{t,x)\'^\v{t,x)\^. 

As in [3] (see also [15] and [25]), we define the followings functional for the 
solution for any t G [0, T*), 



Qo{t) = I / \u{t,x)\ \v{t,y)\ dxdy, 

' x<y 



and 

/oo 
{\u{t,x)\^ + \v{t,x)\^)dx, 
'OO 

Do{t) = / \u{t,x)\'^\v{t,x)\'^dx. 



Qo{t) is called a Bony functional [15] and is similar to the Glimm interaction 
potential. 

Lemma 2.2 There exists constants 6 > such that for the initial data satisfying 
Lo{0) < 6 there holds the following 

^^ + Doit) < 2m(Lo(0))2 (2.5) 

/or t G [0, T^,). Therefore, 

Qoit) + r DoiT)dT < 2miLo{0)ft + QoiO) < 2m(Lo(0))'t + (LoiO))' (2.6) 
Jo 

forte [o,n). 



Proof. By (2.4), we have 
dQoit) 



< -// {\u{t,x)\%\v{t,y)\'dxdy+ \u{t,x)\\\v{t,y)Wdxdy 

J Jx<y J J x<y 

+ {ro{t,x)\v{t,y)\^ + \u{t,x)\^ro{t,y))dxdy 

J J x<y 

oo 



dt 



< -2/ \u{t,x)\^\v{t,x)\^dx 

/oo /"OO /"OO /"OO 

ro{t,x)dx / |vp(i|/+ / |wp(ix / rQ{t,y)dy 
-oo J —oo J —oo J—oo 

/oo 
|w(x,t)|Xx,t)p(ix + 2m(Lo(t))2 
-oo 

= (-2 + 2Lo(0)c)Z}o(t) + 2m(Lo(0))2, 

where we also use Lemma 2.1 in getting last equality. Now, choose a constant 
(5 > such that 

-2 + 2(5c< -1. 

Then, we can get the desired result. The proof is complete. □ 

Next we shall use the above estimates to get the L°° bound of the solution for 
< t < T,. Denote 

Vl{xq, to) = {{x,t) : <t < to, Xo - (to - t) < x < Xq + (to - t)}. 



and 



rij(xo, Xq) = {(x, t) : < t < to, X = Xo + (to - t)}, 
rL(xo, xo) = {(x, t) : < t < to, X = xo - (to - t)}. 



see Fig. 1. 



^L{xo,tQ 



{xo - to, 0) 



(xo, to 




(a;o + to,0) 






Figure 1: Domain Q{xQ,tQ 



Lemma 2.3 For any to G (0,T*), there hold 



I |M|'<g(to), 

JVBixo,ti)) 



rL(xo,to) 



where q{to) = to(mLo(0) + 4cm{Lo{0)f) + 4c(Lo(0))2 + Lo(0). 



Proof. Integrating (2.2) over r2(xo,to) yields 



{\u\^)t + i\u\'^)Jdxdt\= / 2\^{u{miv - iNi))\dxdt 

n{xo,to) ^ ^ "' Jn{xQ,tQ) 

< f f {m{\u\'^ + \v\'^) + 2c\u\'^\v\'^)dxdt 

J jf2(a;o,to) 

< mtoLo(O) + 2c / " D(r)rfr 
7o 

< to(mLo(0) + 4cm(Lo(0))2) + 4c(Lo(0))2 



and 



MvW-i\v\^)x)dxdt\ = 2\^{v{miu - iN2))\ 

fl{xo,to) 



< mtoLoiO) + 2c / " D{r)dr 
Jo 



/O 

< to(mLo(0) + 4cm(Lo(0))2) + 4c(Lo(0))2, 



where we use (2.3) and use Lemmas 2.1 and 2.2. Then the desired estimates 
follow by applying Green formula to the above. The proof is complete. □ 

Lemma 2.4 Fort G [0,T*), there holds the following, 



sup (|n(x,t)| +\v{x,t)\ ) < ( sup(|no(x)| +|i'o(x)| )+2mg(Tj] exp (mT*+2cg(T^)), 
•where q{t) is given in Lemma 2.3. 



Proof. Using the characteristic method, by (2.4) we have 
d 



dt 

and 



|u(t,Xo + (t-to))P < ((m + 2c|vn|wp + m|t;p) 



{t,Xo + {t~to)) 



(t,xa-(t-tn)y 



— \v{t, xo-{t- to))P < ((m + 2c\u\^)\v\'^ + m\u\'^) 
for any xq G R^, to G [0,T^,). Then, 



\u{to,xo)\^ < {\uo{xo -to)|^ + m / \v{t,xo + {r - to))\^dr)hi{to), 

Jo 

\v{to,Xo)\'^ < {\voixo + to)\^ + m / \u{t,Xo - (r - to))prfr)/i2(to), 

Jo 
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where 



/ii(to) =exp(mto + 2c / \v{xo + {t - to),r)\ dr), 

Jo 

/•to 

^2(^0) = exp(mto + 2c / \u{xo - (r - to),r)| dr). 

Jo 

Therefore, these lead to the desired result by Lemma 2.3. The proof is complete. 
D 

Proof of Theorem 1.1. It suffices to get the H^ bound of the solution on 

[0,T*). We differentiate Eqs. 1.1 with respect to x; then 

iiutx + u^x) = -mv^ + {Ni{u,v))^, 

(2.7) 

i{Vtx - Va:x) = -mU^ + (A^2(m, v)),^. 

For the nonlinear terms we have the following estimates by Lemma 2.4: 

for t G [0,T*). Here Ci is a constant depending on m and aj, (3j, 1 < j < 5; 
qi(T^) = (1 + g(T*))^ exp(2mT* + 2cq(T^)). This estimates, together with the 
energy method, we can get the desired bound as (2.1). Thus we can extend the 
solution across t = T^,. The proof is complete. □ 

3 Convergence in L^ 

First we establish the estimates on the difference of smooth solutions. Let («', v') 
be the global smooth solution to (1.1) taking {uq, Vq) G C^{R^) as its initial data 

withJZ.{\u',{x)\'+\v',{x)\^)dx<6. 
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Let U = u — u' and V = v — v'. Then 



Ut + U^ = imV - i(Ni{u, v) - N^^u', v')) , 

^ ^ (3.1) 

Vt-V^ = imU - i{N2{u,v) - N2{u\v')), 

which leads to 

{\U?)t + {\U?)x = m{imVU-i(Nr{u,v) - N^{u\v'))U}, 



(3.2) 
i\y?)t - (l^n. = m{imUV - i[N2{u,v) - iVsK, v'))v} 

and 

m' + \y?)t + {\u?-\y\% = -R (3.3) 

with 

R = ^{2i[Ni{u,v) - Niiu' ,v'))u + 2i[N2{u,v) - N2{u\v'))v}. 
For the nonhnear terms in righthand side in the above, we have the following. 
Lemma 3.1 There exists a constant c* > such that 

\m{imVU - i(Ni{u, v) - iVi(M', v'))U}\ < ri{t, x), 

\^2{imUV - i[N2{u, v) - A^sK, v'))v}\ < n{t, x), 

and 

\R\ < c*r2(t,x,a;), 

where ri(t, x) = m{\U{t, x)p + \V{t, x)p) + c*r2(t, x, x), and 

r2{t,x,y) = \U{t,x)\'{\v{t,y)\' + \v\t,y)\') + {\u{t,x)\' + \u\t,x)\')\V{t,y)\'. 
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Proof. Assumption (Al) implies that there exists a constant c* > such that 



Ni{u,v)-N^{u',v'))U\ = \Ni{u,v)~N,{u',v')\\U\ 

< ^{\u - w'llwp + \u'\\v - v'\\v\ + \u'\\v'\\v - v'\)\U\ 
8 

< -^r2{t,x,x), 



and 



I {N2iu, v) - miu', v'))v\ < -^r^it, X, x) 



4 
Thus, the above estimates give the desired results. The proof is complete. □ 

Now it follows from (3.2) that 

m')t + {\U\%<n{t,x), 

(3.4) 

i\vn-m\<riit,x). 

And to deal with the nonlinear terms in (3.1), we define the following for any 
t> 0, 

Liit)= / i\Uit,x)\'' + \Vit,x)\'')dx, 



<5i(0 = / / r2{t,x,y)dxdy, 

J J x<y 

-^i(^) = / r2{t,x,x)dx. 



Let 



r'o{t,x) = m{\u{t,x)\'^ + \v'{t,x)\'^) + c\u'{t,x)\'^\v'{t,x)\'' 



and 



Q',{t)= / / \u'{t,x)\'\v'{t,y)\'dxdy, 

J Jx<y 
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and 



oo 



Lo(t) = / {\u'{t,x)\'^+\v'{t,x)\'^)dx, 

J— oo 

roo 

Do{t)= \u{t,x)\'^\v'{t,x)\'^dx. 



Lemma 3.2 There exist constants 6 > and K > such that if Lq{0) < 6 and 
L()(0) < 6 then 

j^{L^{t) + KQ^{t)) + D^{t) < (2mLo(0)+2mL[,(0) + cDo(t) + cZ};,(t))Li(t) (3.5) 

for t > 0. Therefore 

Li{t) + KQi{t) < (Li(0) + KQi(O)) exp(/i3(t)) (3.6) 

and 

I Di{r)dT< {Li{0) + KQi{0)){l + {Am6 + S^ + 2mSH) f exp(/i3(r))rfr) (3.7) 
Jo Jo 

for for t > 0, where 

hsit) = 2mLo(0)t + 2mL'{0)t + f {cDo{r) + cD'{T))dT. 

Jo 

Proof. (3.4) yields that 

^Li(t)<2c,Di(t), 
dt 

and 

d 



-rMt) < -2D,{t)+ r,{t,x){\v{t,y)\' + \v'{t,y)\')dxdy 

at J Jx<y 
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+ / / \Uit,x)Wro{t,y) + r'Q{t,y))dxdy 



x<y 



{\u{t,x)\ +\u{t,x)\ )ri{t,x)dxdy 

x<y 



{ro{t,x) + r'Q{t,x))\V{t,y)\ dxdy 

lx<y 

< -2Di(t) + (mLi(t) + c,Di(t))(Lo(t) + L[,(t)) 
+ (mLo(t) + mLo(t) + cD^it) + cD[,(t))Li(t) 

< (-2 + c,(Lo(0) + L[,(0)))Di(t) 

+ (2mLo(0) + 2mL[,(0) + cDo(t) + cD'^it)) Li{t) , 

where Lemma 2.1 is used in the last inequahty. These estimates give (3.5) for 
small Lo(0) < 5 and -^q(O) < 5 and for Large K, where 5 = ^. (3.6) is a 
consequence of (3.5). 

To prove (3.7), we integrate (3.5) over [0,t], then 

|*Di(r)dr < Li(0) + Kgi(0) + y'*(4m5 + cDo(r)+cD[,(r))Li(r)rfT 

< (Li(0) + KQi{0)) + {4m6 + (5^ + 2m6H) f Li{T)dT 

Jo 

< {U{Q) + KQi{Q)){l + {Am6 + 5^ + 2m5H) f\xp{h3{T))dT), 

Jo 

where Lemma 2.2 and (3.6) are used. The proof is complete. □ 
Lemma 3.3 Let h^i^t) be given by Lemma 3.2. Then for t > 0, 

h^{t) < 2m(Lo(0) + L[,(0) + {L,{Q)f + {L',{Q)f)t + {L,{Q)f + (L[,(0))^ 
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Proposition 3.1 Let {u,v) and {u',v') be two classical solutions to 1.1 with ini- 
tial data {uq,Vq) and{uQ,v') respectively, and suppose that \\uo\W2rjii) + \\vo\\'j^2rjii\ < 
d and ||wol|£,2Cj:ji-) + ||'L'q||^2(^i-) < S. Then, there exist constants, C\, C2 and Cj,, de- 
pending on the ^, such that for t >0, 

\\u{t) -U'{t)\\l2^jii) + \\v{t) -V'{t)\\l2^jii) < (||Wo-<||i2(^i) + ||t;o-^^ollL2(i?i))^4(t), 

where h^i^t) = c\ exp(c2t + C3) . 
Proof. Indeed, we have 

Qi(0)< / (|f/(0,x)|2+|\/(0,x)|2)dx(Lo(0) + L'o(0)), 



which, together with (3.6), gives the proof of this proposition. Thus the proof is 
complete. □ 

Proof of Theorem 1.2 By Proposition 3.1, we have 



\Uk{t) - Mj(t) 1 112(^1) + \\Vk{t) - t'i(t)||i2(^i 



< 



UkO — ^io| 1^2(^1) + ll^fcO — Vjo\\L2(Ri))h4(t), 



which imphes that there exists a unique (mooi "^oo) ^ ^{[0, +00), L^R^)) such that 

ll^fc — '"oo||c([0,r],L(R2)) + Ikfc — 'i'oo||c([0,T],L(i?2)) ~^ 

for any T > 0. 

Now to prove that (Mooj'^oo) is a weak solution of (1.1), let 

Uk,j =Uk- Uj, Vk,j = Vk- Vj, 
16 



and 



rk,j = \Uk,jit,x)\'^{\vkit,x)\'^ + \vj{t,x)\'^) + {\ukit,x)\'^ + \uj{t,x)\'^)\Vkjit,x)\'^. 
By (3.7), for any T > 0, there exists a constant Ci(T) > such that for / = 1, 2 



I Ni {uk ,Vk) - Ni {uj ,Vj)\dx < c^ Tkjdx 

< C,iT)cM\Uk,mlHm) + llVkAmUm)) 

< C^iT)c4\\Uk,M\\hiRr) + ll^.,,(0)||i.(^.)), 

which leads to the strong convergence of the nonhnear terms. Therefore, (uoo, "^oo) 
is a weak solution of (1.1). The proof of Theorem 1.2 is complete. □ 
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